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Abstract* Previously, the stability of a premixed flame in a stagnation flow was discussed for a res- 
tricted class of disturbances that are self-similar to the basic undisturbed flow; thus, flame fronts with corru- 
gations only in the cross stream direction were considered. Here, we consider a more general class of three- 
dimensional flame front perturbations which also permits corrugations in the streamwise direction. It is 
shown that, because of the stretch experienced by the flame, the hydrodynamic instability is limited only to 
disturbances of short wavelength. If in addition diffusion effects have a stabilizing influence, as would be 
the case for mixtures with Lewis number greater than one, a stretched flame could be absolutely stable. 
Instabilities occur when the Lewis number is below some critical value less than one. Neutral stability boun- 
daries are presented in terms of the Lewis number, the strain rate and the appropriate wavenumbers. Beyond 
the stability threshold the two-dimensional self-similar modes always grow first. However, if disturbances of 
long wavelength are excluded, it is possible for the three-dimensional modes to be the least stable one. 
Accordingly, the pattern that will be observed on the flame front, at the onset of instability, will consist of 
either ridges in the direction of stretch or the more common three-dimensional cellular structure. 
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INT H O PWTXON 

The response of a premixed flame to stretch is a problem of fundamental 
importance in combustion theory. Results of studies on the structure, 
propagation and extinction of stretched flames have been used in experimental 
and theoretical investigations (Wu and Law 1984; Law et al. 1986; Tien and 
Matalon 1990; Dixon-Lewis 1991) aimed at understanding the structure and 
burning characteristics of laminar flames. Results of these studies have been 
also applied to turbulent combustion, in particular, in the "reaction sheet" 
regime (Williams, 1985; Peters, 1986) where the turbulent flame is regarded as 
an ensemble of stretched laminar flamelets. Thus, considerable attention has 
been given to stretched flames in the literature where a large number of 
theoretical and experimental papers can be found. In spite of these numerous 
investigations the effect of stretch on flame front instabilities has not been 
completely clarified. 

Experimental studies (Ishizuka and Law, 1982; Ishizuka, 1988) on the 
stability limits of premixed flames clearly indicate the importance of flame 
stretch. Observations were made on propane/air premixed flames in an 
axisymmetric stagnation point flow. The flow field in such a configuration 
imposes on the bulk of the flame a well defined strain rate, or stretch, 
controlled by the mass flow rate of the incoming fresh mixture. Depending on 
the strain rate, various patterns appear on the flame front for rich 
mixtures. For a weak strain rate a three dimensional cellular flame evolves, 
as may be expected; cellular flames are commonly observed in experiments 
performed in large tubes (Clavin, 1985), namely in the absence of stretch. 

Upon increasing the strain rate, the flame approaches the stagnation point and 
a different pattern, characterized by ridges along the radial direction, 
emerges on the surface of the flame front. For sufficiently large values of 


the strain rate, all instabilities are suppressed resulting with a smooth flat 
flame surface. In contrast, for lean mixtures the observations indicate an 
absolutely stable flame front for equivalent values of the strain rate. These 
experiments clearly point out the role of flame stretch, combined with 
hydrodynamical and diffusional- thermal effects, in suppressing/promoting 
instabilities . 

A complete theoretical study on the effect of stretch on flame stability 
has not yet been given. Although there Is no difficulty in formulating the 
linearized disturbance equations for a flame in a stagnation-point flow, the 
resulting problem consists of a coupled system of partial differential 
equations with variable coefficients. These coefficients, which depend on the 
unperturbed base flow, can be written explicitly only for a hydrodynamic 
model, where the flame is treated as a surface of discontinuity (Eteng et al., 
1986). Furthermore, since the basic state is a non-parallel flow, with 
velocities varying normal and tangential to the flame front, the disturbance 
equations are not separable. Thus, the linear stability problem cannot be 
resolved using normal modes in both directions that are parallel to the flame 
front. An exception is a class of disturbances self -similar to the basic 
state which has been considered by Kim and Matalon (1990). In this case, 
perturbations have periodic waviness in the cross -stream direction only with 
no corrugations in the direction of flow divergence. Self similar 
disturbances were first considered by Gortler (1955) and more recently by 
Wilson and Gladwell (1978) , In analyzing the stability of the incompressible 
viscous stagnation-point flow without chemical reaction. More general 
disturbances were later admitted by Spalart (1989) and Brattkus and Davis 
(1990) in dealing with the fluid dynamical problem. Spalart solved the 
initial value problem numerically and found that the long term behavior of the 
solution approaches the spatial structure of the self similar mode. Brattkus 
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and Davis expressed the disturbances in terms of an infinite series of coupled 
modes and concluded from their analysis that the least stable mode is that 
corresponding to a self- similar disturbance. Although these studies provide 
some justification for the use of self-similar disturbances, perturbations 
with periodic waviness in the streamwise direction are not of negligible 
significance, especially for the combustion problem. These modes are 
undoubtedly important in the nonlinear development of flame disturbances and 
must be included in any description of the three dimensional cellular pattern 
which has been observed in experiments. 

In addition to the analysis of Kim and Matalon, the only other reported 
study on the stability of a flame in a stagnation-point flow is that of 
Sivashinsky et al . (1982). Their analysis, however, is restricted to a very 
weak stretch rate. Furthermore, they suppressed the hydrodynamic instability 
by neglecting the effect of thermal expansion, thus deriving results that 
pertain to the effect of stretch on the diffusional- thermal instability only. 

In the present work, we analyze the stability of a premixed flame in a 
plane stagnation-point and consider a more general class of three-dimensional 
disturbances which, in particular, allow for flame front corrugations in the 
streamwise direction. As is well known, the hydrodynamic instability (Landau, 
1944; Darrieus, 1945) does not permit the existence of premixed flames that 
are too flat. Although diffusion often introduces stabilizing influences, as 
is the case for mixtures with Lewis number greater then one, long wavelength 
disturbances are hydrodynamically unstable. Our present aim Is to determine 
whether a sufficiently large stretch can stabilize these long wavelength 
disturbances. We have also Incorporated thermo-diffusive effects in our model 
in order to determine the extent that a flame needs to be stretched to remain 
stable. For two-dimensional self -similar disturbances, Kim and Matalon (1990) 
found that stretch indeed stabilizes the flame. In particular, they provided 


the modification resulting from stretch, to the neutral stability curve drawn 
previously for freely propagating flames (Matalon and Matkowsky, 1984; Jackson 
and Kapila, 1984), thus identifying the critical conditions at the stability 
threshold. Here, we shall not only extend their results, but also identify 
the effect that such disturbances have on the pattern that evolves on the 
flame front and in the induced flow field when an instability develops. It is 
shown that beyond the stability threshold, the two-dimensional self-similar 
modes always grow first. However, if disturbances of long wavelength are 

i 

excluded, it is possible for the three-dimensional modes to be the least 
stable one. Accordingly, the pattern that will be observed on the flame front 
at the onset of the instability, will consist of either ridges in the 
direction of stretch, or the more common three-dimensional cellular structure. 
These results appear in agreement with the experimental observations 
summarized above, though a direct comparison requires analyzing the 
axisymmetric configuration. Having established a procedure for testing self- 
similar base flows for stability against general disturbances, we shall 
discuss the stability of a premixed flame in an axisymmetric stagnation-point 
flow in a future publication. 
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FORMULATION 

The mathematical formulation is identical to that of Kim and Matalon 
(1990). A combustible mixture impinges toward a flat surface held at x - 0 
(see Fig. 1). Far upstream, the flow field is the classical stagnation-point 
flow characterized by the strain rate e. The flame, which separates the 
burned products from the fresh unburned mixture, is considered thin and is 
therefore represented by the surface $(x,y,z,t) — 0, where ^ > 0 is the burned 
gas region. The flame structure resolved on the scale of the diffusion length 
Lp - D T /S f , where D T is the mixture thermal diffusivity and S f the laminar 

flame speed, yields jump conditions for mass and momentum across the flame, as 
well as an equation describing the evolution of the flame front (Matalon and 
Matkowsky, 1983). Since viscosity has only a secondary effect on flame 
stability, it will be neglected everywhere. Consequently, the flow field on 
either side of the flame is described to all orders by the inviscid, 
incompressible Euler equations but with distinct densities in the unburned and 
burned gas. The jump conditions and the flame speed equation are only correct 
to 0(5), where 5 - I^/L « 1 is the ratio of the diffusion length to a 

hydrodynamic length L. 

In writing the governing equations the length L is used as a unit of 
distance, the laminar flame speed as a unit of speed, L/S^ as a unit of 

time and the state of the fresh mixture p Q , p Q , T Q is used to 

nondimens ionalize the density p, pressure p and temperature T, respectively. 
These equations are: 
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for $(x,y,z , t) < 0 


V*v - 0 


across tf(x,y,z,t) - 0 

[p (v*n) - V f ] - 5 (£ n a ) K (3) 

f 

[v x n] - - 6/9-j^Vxv] + 2(<7-l) VxnJ’ (4) 

[p + p (v»n) (v*n - V f )] - - (ct- 1)(1+0) V»n + V f (in a) kJ- (5) 

where [•] denotes the jump in the quantity, i.e. (•). .'-(•) 

J * * v 'burned v ''unburned’ 


at tf(x,y,z,t) = O' 


v*n - V f - 1 - 5(/3+ 7 i/2) K 


In the above D/Dt is the convective derivative and the variable p stands for 
the small pressure deviations from the ambient pressure level (of the order of 
a squared Mach number) . The density is given by 


r 1 * < 

* 

1/a * > 


where a is the thermal expansion parameter (a > 1). The velocity of the flame 
surface and the unit normal vector n are 

V _ . _L_ 21 - _ _v*_ / q \ 

f |v*| at ’ |v*| • 

The flame stretch is given by 


K = V«n - n* Vx(vxn) 


(9) 


KIM Hill i INI II llll IIMI Nil II I mill II IN Hi 
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with the right hand side evaluated at tf(x,y,z,t) - O'. The remaining 
coefficients are 

0 - , 7 - d€. i - (E/RT q ) (Le - 1 )/ct 2 

where E is the activation energy of the chemical reaction and R the gas 
constant. Note that /3 and 7 are both positive quantities which depend solely 
on the thermal expansion parameter, and I is a scaled Lewis number that 
measures the deviation of Le from one. Finally the boundary conditions are 

u = 0 at x = 0 (10) 

v - ( - ex , ey , 0) as x -*■ -«> (11) 

where v - (u,v,w) and e - (L/S f )e is the dimensionless strain rate. 

The basic state corresponds to a steady planar flame front located at 
x - -d. The velocity field is given by (Eteng et al., 1986 ) 

v = ( F(x), -yF'(x), 0 }, ( 12 ) 


where F(x) - F Q (x) + 5 F^(x) and prime denotes differentiation with respect to 


x. Here, 


F q (x) 


F 1 (x) 


- e (x + a Q ) -*> < x < -d 

sin qx -d < x < 0 

o-l ^ 


- ea, 


„ ol/_2_vl/2 
£a 2 ( a-l> 


sin qx - qx cos qx 


-00 < x < - d 
-d < x < 0 


where q » ejo-l/o , The pressure field takes the form 

- 2^0 + € V) - 5 e 2 a 1 (x + a Q ) -<*> < x < d 

p - « ( 13 ) 

; k[ F o + «V] - % F i - -d < x < 0. 
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j 




The flame standoff distance d and displacement distance a are expressed as 
d - dp + Sd-^ and a - a^ + S a^, where 

d q - rj/c , d 1 — - - y^rf/2 

a o “ d o ’ 1//€ 9 a l “ d l + ^ + 7 ^/ 2 

with rj - (cr tan 1)/Ja- 1 . Note that for a very weak stretch (c « 1) the 
flame recedes back to the nearly uniform flow, far upstream. For 0(1) values 
of € the flame stands at a finite distance from the plate, thus justifying the 
inviscid approximation assumed earlier. Only when c » 1 does viscous effects 
become important because the flame, in this limit, is pushed into the thin 
viscous boundary layer adjacent to the plate. 



in mi min I'm rn min ii mi imp iii'iiiniii ill iii'iin 
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r U(x,y)’ 

V 


v(x,y) 

w 


w(x,y) 

. p . 




DISTURBANCE EQUATIONS 

We introduce small disturbances, superimposed to the basic state (12)- 
(13), of the form 


, ikz+wt (14) 


and express the perturbed flame front as 

x — -d + #( y)e^ Z +° )t (15) 

where k is the wavenumber associated with the spatial periodicity in the z 
direction and to is the growth rate. Note, that using (14) -(15), a solution 
that assumes normal modes in z and t, but not in y, has been sought. 
Substituting into the governing equations (l)-(2) and then linearizing about 
the basic state, one obtains 


Fv 


x 


u + v + ikw - 0 
x y 


Fu x - yF'u y + (co + F')u - - p x /p 

yF'v y + (» - F')v - yF»a - - p y /p 

Fw - yF'w + tow - - ikp /p 
x y 


(16a) 

(16b) 

(16c) 

(16d) 


where, here and thereafter, subscripts denote partial differentiation. 

In the absence of vorticity far upstream, and as a result of the density 
being constant, the flow field up to the flame can be determined from a 

2 

velocity potential <p(x,y,x,t) which satisfies V <p - 0. Writing 

a , x ikz+tot 
<p * 9 ?(x,y) e 

Laplace's equation yields 


A , A i A 

<p + a? - k <p 
xx yy 


( 17 ) 
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from which the velocity field can be determined as 

’ v “ . V - ik* . (18) 

The pressure is obtained from any one of the components of the momentum 
equation, giving 

P - yF '* y - F* x - W*. (19) 

Thus, disturbances in the unburned gas are determined in terms of a single 
variable cp. In the burned gas, however, the flow, is no longer potential 
because the vorticity produced at the flame is transported downstream. 
Consequently the stability problem reduces to solving (17) in the unburned 
region and the system (16) in the burned region. 

Expressing the jump conditions (3) -(5) and the flame speed relation (6) 
at the mean position x - -dg, then linearizing about the basic state and using 

the solution expressed in terms of <p for the unburned region, one obtains the 
following conditions: 


u = * x - «(a-l)|^+(l+» ? )^j 1$^- k 2 $j + ^e 2 »;( y#) y + k ^]} ( 20a ) 

v - *y - 2^»?|j“(l+^/»7)y0 + [u+e-eo/ri)$ y 

+ e y$yy + ( w+e )^y + <y* y } (20b) 

w = ik* - ik(o-l)$ - ik ^^(w- ta/r\')^> + ey$ y + w* + ey*yj- (20c) 

P “ -«y£ y - l-«e(l-i7)^j* x - oxp + 6(o-l)^l+)9+^(»7+2)J - k 2 ^ 7 fj 

+ 7-f ^yy - k 2 *j + ^j^(2e 2 -w 2 )^ - 2ew y# - f2 y 2 ^y y ]} ( 20d ) 
V x “ («*+«)$ + <y$ y + 5 {’^ + 2 ^)[^yy ‘ fc2 #] + 2 ^(^'l)[^yy - k ^]} (20e) 


(20c) 


+ ^(t-D^yy * k ^]} 


(20e) 


» m Him i hi ilium mill i jii< illinium iinuiiiHiiiiiiiiijiiWiiiiiA iwm min iiini m mi iimm 
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Far upstream, i.e. as x -+ -°o, all velocities must vanish and therefore $> must 
be bounded. At the plate, only the axial velocity must vanish, i.e. u - 0 at 
x - 0, as appropriate for inviscid flows. 

The system (16) for x > -d together with equation (17) for x < -d, the 
relations (20a-e) at the flame front and the boundary conditions at x - •« and 
at x - 0, constitute an eigenvalue problem for the determination of the growth 
rate The eigenvalues u> and the associated eigenfunctions depend on the 
following parameters: the thermal expansion a, the strain rate e, the 
(reduced) Lewis number i, the flame thickness 8 and the form of the 
disturbance represented by k, the wavenumber associated with the z- direction, 
and an appropriate parameter that remains to be determined, identifying the 
form of the corrugations in the y-direction. 


FORM OF THE DISTURBANCES 


As already pointed out in the introduction, the eigenvalue problem we 
have just obtained is not separable so that normal modes in y are not 
admissible solutions. An exception is the self similar form 


{ u, v, w, p } - { U(x) , yV(x), W(x) , P(x) } 


( 21 ) 


with $ - const., considered by Kim and Matalon (1990). In this case, all 
variables are functions of the axial coordinate x, with the exception of the 
streamwise velocity component which is linearly proportional to y. Clearly, 
this follows from the fact that the stagnation-point flow is a local solution 
valid near the axis (y - 0). It is plausible, therefore, to expect that more 
general disturbances will include higher powers of y as well. 

Following Brattkus and Davis (1990), further insight in choosing the 
allowable disturbances can be obtained by examining the vorticity equation 


+ F fl - yF' 0 - - ikyF" w + F" v + yF' " u , 

x y 


( 22 ) 


A 


obtained by cross-differentiating equations (16b) and (16c), where Q - 

is the z- component of the vorticity perturbation. Note that the base flow 
contains vorticity in the burned gas = -yF M , which gives rise to the 

nonzero terms on the right hand side of equation (22) . If one ignores 


A 


u 


y 


and 


considers a potential flow F = -x, the right hand side of equation (22) 
vanishes. The simplified equation, for -d < x < 0, < y < «>, admits 

a n+cj n 

separable solutions of the form ft - x' w y . Depending on n, and on whether 
the disturbances grow or decay, the vorticity may either vanish or amplify as 
x -+ 0; in the latter case the vorticity decays in the boundary layer near the 
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wall due to viscous action. In any case 0 increases with distance from the 
n 

axis as y . 

In the presence of vorticity in the base flow, equation (21) is no longer 
separable and therefore the disturbances are not expressible In a simple form 
as in (22). Nevertheless, we may write 


u - S U (x)y n , 
n=0 


v - S V (x)y n , 
n-0 


w — E ikW (x)y 
n-0 


(23) 


£ - E $ (x)y n , # - 2 A n y n 

n-0 n-0 

expecting that now the various modes, corresponding to different indices n, 
are all coupled. Brattkus and Davis (1990) have also expressed the 
disturbances as an infinite series, but in terms of Hermite polynomials which 


arise naturally when viscous effects are incorporated in the vorticity 
equation. Finally, it should be pointed out that unlike the self-similar 
disturbances (21), which are essentially two dimensional with corrugations 
only in the z- direction, the disturbances considered here are three 
dimensional . 
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Upon substituting (23) into the governing equations (16) - (17) and the 
boundary conditions (20), the following system results: 
for x < -d„: 


*n ' k $ n “ -(n+D(n+2)$ n+2 


for -d^ < x < 0: 


U' + (n+l)V - - k W - 0 
n n+1 n 


FU' + [w - (n-l)F' ]U - - a?' 
n n n 


FV^ + [w - (n+l)F']V n - - <7(n+l)P n+1 + F"U n l 


FW' + [u> - nF' ]W - - ctP 
n n n 


U n - - 6 (cr- l)^/3+(l+?7)^J j^(n+l) (n+2)A n+2 ’ k ^ A n] + 2 ^ £ ^(n+l)A n 

+ ^[(n+l)(n+2)$ n+2 - k 2 $JJ- (26a) 

3 

V n - + € ^ A n-l] ’ ^{^ (!«/■>)*„. i 

+ [w+£(n+l)-eo/rj)](n+l)A n+1 + [«+(n+l) £ ] (n+l)$ n+] J (26b) 

W n “ $ n‘ (a ' 1)A n" 2 ^^[w+en-£a/ij) ]A n + (w+nO^J- (26c) 

P n ” -( w+ne )^ n - l-$£(l-»?)^j*^ + 5(a-l)^7ij^(n+l)(n+2)$ n+2 - k 2 $J 

+ [l+^(r ? +2)j [(n+l)(n+2)A n+2 - k 2 A n ] 

+ ~!^r\ j^2€ 2 -w 2 -2n£Ci>- £ 2 n(n-l) j A^|- (26d) 


P — -(t<;+ne)$ 
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n 


[w+e(n+l)]A n + 5«J- p+^ (n+l)(n+2)A n+2 - k 2 A n j 


+ J 1 ( n -l)[(n+l)(n+2)$ n+2 - k 2 *J)- (26e) 


at x - 0 


U - 0 
n 


(27a) 


atx-1 -°° 


$ are bounded 
n 


(27b) 


We note that the problem for the n order mode is coupled to the systems 
corresponding to the n-1, n+1 and n+2 modes. The mode coupling in the 
equations appears only on the right hand side of equation (25c) and is a 
result of downstream pressure gradients, and of vorticity produced by the base 
flow at the flame front. Mode coupling in the conditions (26) is of two 
origins: The first, which is a result of vorticity generation at the 

perturbed flame front, appears in the 0(1) terms on the right hand side of 
equation (26b) . The second, which is a result of flame front curvature and 
flame stretch (effects that are of the same order as the flame thickness) , 
appears in the 0(5) terms of these relations. 

A direct observation suggests that the eigenfunctions of the system (24)- 


(27) fall into two distinct classes, with the velocity v being either 
symmetric or anti- symmetric about the axis y “ 0. Each of these two classes 
has its own set of eigenvalues. We refer to the even or odd solution as the 


series containing the sequence 


{ 


U 


V, 


2n’ 2n’ 2n-l 


W, 


2n 


n = 0,1,2. 



or, 


{ 


A II V W 

2n-l’ 2n-l’ 2n’ 2n-l 


n - 0,1,2... I*, respectively. The classification 


is made according to whether the disturbed flame front # and hence u and w are 
even or odd functions of y; v always has the opposite parity. We determine 
these solutions by truncating the series at n - N and finding the eigenvalues 
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at any given level of truncation. For a given N, one of the two following 
system of equations results for the determination of the eigenvalues: 


System I 

i 

U N-1 + ' k2 Vl “ 0 (28a) 

FU' i + [« - (N-2)F']U n1 - - aP^ (28b) 

FV^ + [« - (N+1)F']V n - F"U n1 (28 c) 

i 

FWn_i + [w - (N-1)F']W N1 - - <7P n1 (28d) 

System II 

U N - “ ° (29a) 

FU^ + [w - (N-1 )F']U n - - aP^ (29b) 

FV^i + [» - NF']V n1 - - aNP N + F"U n _ 2 (29 c) 

FW ^ + [» - NF']W n - - aP N (29d) 


The remaining equations, for the variables with lower indices, are only needed 
for the determination of the eigenfunctions. We note parenthetically that 
equation (29c) is uncoupled from the remaining equations in system II and is 
therefore not needed for the determination of w. Each of these two systems is 
to be supplemented by the conditions (26a-d), truncated appropriately. The 
latter requires the determination of and however, from equation (24) 

one notes that the solution that decays as x -+ always satisfies 

which, when substituted in (26e) determines 

As N increases from N - 0, we first find the solution - 0; for N - 1 

we find the first even mode (Uq ; ; Wq) which is the self similar solution 


n |! hi ii ii i mu mil ii 1 ii numRiiK iniiimi i I ! ltn f mi llliiiimMipni mu him i iim miiniii imt n nm i nil iiiiihiiiiiiihiiiri mm iiiim hi iniiiKim n iliMUMII 'irill 1 nil I ml he i i i sw rum iummmi i mill 



17 


(see equation 21) and the first odd mode (U^; Vq; W^) ; for N — 2 we find the 

second odd and even modes (U^; V^; W^) and (I^, U^; V^; , ^q) 

respectively, etc,.. . In general, when N Is odd system I determines the 
eigenvalues of the N-even mode and system II determines the eigenvalues for 
the N-odd mode; when N is even system I determines the eigenvalues of the N- 
odd mode and system II determines the eigenvalues for the N-even mode. 
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STABILITY OF A STRETCHED FLAKE 


| 


The mathematical formulation presented above is valid for 8 « 1, so that 
the solution must be expressed in power series of 6 . Here we only follow the 
first two terms in this expansion and, in particular, express the dispersion 
relation as 

co « coQ(a,€;k,N) + 6co^ (a , e , i ; k , N) . (30) 

For a given N, the eigenvalues for the even and odd modes are obtained by 
solving equations (28) or (29) with their appropriate boundary conditions. 

The leading problem, obtained by setting 8 * 0, yields co^. The next term co^ 


is obtained as a solvability condition from a system with the same 
differential operator as that of the leading problem but with nonhomogeneous 
terms that depend on the leading order solution. The solutions presented 
below were obtained by numerically solving these systems using the boundary 
value solver COLSYS (Ascher et al . , 1981). 

The dependence of the solution on the strain rate e can be scaled out in 
the whole problem, as shown by Kim and Matalon (1990). Consequently, the 
dispersion relation (30) can be rewritten as 


Cd 


eaQ(a;/c,N) + 6e (a , i ; * , N) 


(31) 


where k, - k/e . 

The leading order term in the dispersion relation corresponds to the 
modification of the Darrieus -Landau instability resulting from flame stretch. 
For very small e, the flame is located far upstream in a nearly uniform flow. 
Thus, the hydrodynamic disturbances grow in the infinitely long downstream 
region of burned gas according to 
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(32) 
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Landau, 1944; Darrieus , 1945). For moderate values of e, the flame Is located 
at the mean position x — -d. Consequently the hydrodynamic disturbances which 
are now confined to the region -d < x < 0 of burned gas, must vanish at a 
finite distance from the flame front. This damping effect is one source of 
stabilization, the only one experienced by the two dimensional self- similar 
disturbances. An additional source of stabilization comes from the stretching 
that three dimensional disturbances experience in the direction of flow 
divergence. These two effects modify the growth rate (32), as shown in 
Fig. 2, where the leading order term of the dispersion relation is plotted 
against k for a selected number of even and odd modes corresponding to various 
N. [Here, and in all following figures, solid curves result from solutions of 
system I while dashed curves result from solutions of system II; the notation 
associated with the different modes consists of a digit which stands for the 

value N and the letters E for even and 0 for odd] . We note from Fig. 2a that 

disturbances with long wavelength in the z- direction (i.e. with small k) are 
now damped out. However, for a given k, disturbances with short wavelength in 
the y-direction (i.e. corresponding to large N) are the more stable one. For 

a fixed e there exists therefore a value k" (N) such that the only growing 

modes are those with k > k (N) , with k increasing with N. For comparison, we 
have plotted in Fig. 2b the relative growth rate o>q compared to that of an 

unstretched flame cjq^ * For a fixed k this graph identifies the dependence of 

the growth rate on 1/e . Note that the ratio is less then one for 

e > 0 and approaches one as e -+ 0. 

The effects of diffusion appear in the 0(6) term of the dispersion 
relation (30). For an unstretched flame, the correction to the growth rate 


is of the form 


.( 0 ) 


f ri g(g-l) (l+£) + 2a/?(l+a 0 )- | J-l 

2 ct+(ct+1)q 0 J + L 2 
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a+(<7+l)a Q 


(33) 


(Pelce and Clavin, 1982; Matalon and Matkowsky, 1983). In the presence of 
stretch it is clear from the boundary conditions (26) that the solution also 
depends linearly on £ so that one can write 

— q^(<7;k,N) + ia^2(g>' c >^) • (32) 


Our calculations indicate that, similar to the case e — 0, the coefficients 
a u and a 12 are both negative (see Fig. 3). This implies that diffusion has 

always a stabilizing influence when the Lewis number Le is greater than one 
(i.e. £ > 0) . Diffusion is a source of instability only when Le < 1; the 
precise determination of the critical £ will be discussed later. Note that 
the coefficients and c* 12 are quite insensitive to variations in N; for N 

as large as 5 the variations in remain indistinguishable while those in 


a^ 2 are less than 10%. 

The combined effects of hydrodynamic, including stretch, and diffusion 
can be understood by considering the expansion (30) as an approximation to the 
exact dispersion relation. This procedure was found to provide qualitatively 
correct results for the case e - 0 (Matalon and Matkowsky, 1984) which were 
also quantitatively in agreement with the numerical results of Jackson and 
Kapila (1984) for small k. The results for a stretched flame are shown in 
Fig. 4 where the growth rate w is plotted versus k for a selected number of 
modes N. The stabilization effect of diffusion is clearly seen when compared 
to Fig. 2. First, disturbances with short wavelength in the z-direction are 

damped out so that, for a given e, growing modes are limited to k < k < k . 

Second, both k' and k + decrease with increasing N so that, for sufficiently 
large N, all modes are damped out (i.e. w < 0 for all k and N) resulting with 
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a stable flame front. We note from Fig. 4a, which corresponds to £ - - 0.53, 
that the most dangerous disturbance is IE, namely the two-dimensional self- 
similar disturbance, with 0.5 < k < 3. However, for £ — - 0.54 (see Fig. 4b), 
three dimensional disturbances with N < 2 and relatively much larger k, grow 
at a faster rate than mode IE does. As £ decreases further, more of the 
three-dimensional disturbances grow faster than IE; for example, when 
£ - - 0.6 all five modes shown in Fig. 4b have growth rates that are larger 
than that of IE. The implication of this observation will be discussed next. 

The neutral stability boundary is obtained by setting w - 0 in the 
dispersion relation (31). The results, plotted in Fig. 5, are presented in 
the k-2 plane for several values of the parameter eS. The various curves 
correspond to different values of N. The region to the right of these curves 
is where w is negative and the flame is stable. The stability boundary is the 
envelope of these curves and the critical £ c is defined as the maximum value 

of £ taken over all /c and N. Note that £ < 0 so that a stretched flame in a 

c 

premixed combustible mixture with Lewis number Le > 1 , such as a lean 
hydrocarbon/air mixture, is stable. Upon decreasing i, for example by making 
the mixture richer, an instability develops. From Fig. 5 one sees that at the 
stability threshold the two-dimensional self-similar disturbances grow first 
and that the transition occurs at a larger value of i, the smaller eS is. The 
pattern that will be observed on the flame front, in this case, is one with 
ridges in the z- direction only. If the finite width of the apparatus in the 
laboratory exclude disturbances with long wavelength, it is possible for a 
three dimensional mode to be the least stable one and, consequently, for a 
cellular flame to be observed at the onset of instability. For example, If 
disturbances with k < 3 and N < N° are excluded for eS - 2, one of the dashed 

^ ■ji' 

curves characterized by some k and N > N° will be crossed first upon 
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decreasing i. Thus, depending on the strain rate, two or three dimensional 

patterns could result on the flame front of rich hydrocarbon/air mixtures . I 

- • | 

The structure of the flow field at the onset of instability is quite different 

in these two cases, as shown in Fig. 6, In the figure, streamlines of the 

perturbed flow field on either side of the flame are plotted in the y plane. I 

Case (a) corresponds to the IE solution derived from system I (eqs. 28) 

associated with the right most point of the solid curve in Fig. 5; case (b) 1 

corresponds to the 10 solution derived from system II (eqs. 29) associated 

with the right most point of the dashed curve in Fig. 5. We note that in the 

i 

latter case the flow pattern consists of several arrays of vortices stacked in | 

the upstream direction and periodically distributed in the spanwise direction 
z. This whole structure is convected along the diverging streamlines of the 
mean flow. 
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CONCLUSIONS 


We have considered the stability of a premixed flame in a stagnation- 
point flow against a class of three dimensional flame front perturbations. It 
has been concluded that the least stable mode corresponds to disturbances 
which are self -similar to the basic state, in accord with the known stability 
results for a viscous incompressible stagnation-point flow without chemical 
reactions. If disturbances of long wavelength are excluded, however, it is 
possible for a three- dimens ioanl mode to grow first. This of course 
determines the pattern that is likely to be observed on the flame front at the 
onset of instability. We have also determine in this study the stability 
criteria in terms of the strain rate, the Lewis number and the thermal 
expansion parameter. 
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Streamlines of the perturbed flow field, on either side of the flame 
located at x/d - -1, at the onset of instability; (a) from the 
solution IE derived from system I for a - 5, e6 “0.8, i - - 0.43 
and k - 0.8; (b) from the solution 10 derived from system II for 

a — 5. e 8 — 0.8. i — - 0.53 and k — 4.4. 
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